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1. The WHIM 

Baryonic matter - protons, neutrons, and electrons - belongs to a class of massive elementary 
particles which make up the atoms of all materials on Earth and in the stars. From the Big Bang 
theory and the results of WMAP, baryonic matter contributes no more than 5 percent of the cosmos, 
though it remains not completely accounted for. From a census of the observable baryonic density 
in the current Universe (e.g. Shull, Smith & Danforth (2012)), all the stars, dust and gas within 
galaxies constitute less than 50 percent of baryons predicted by the theory. By analyzing light from 
distant quasars, the number of baryons in cosmic history can be counted. Its amount has remained 
constant even going back to 10 billion years ago. Quantifying the missing baryons is essential to 
obtaining a complete picture of galaxy formation and evolution. The missing baryons are generally 
thought to reside between galaxies, as warm-hot intergalactic medium (WHIM) with temperatures 
of 10 5 — 10 7 K. Material in the WHIM is highly ionized and can only be observed at far-ultraviolet 
or low-energy X-ray wavelengths. This diffuse plasma is a common prediction in the A cold dark 
matter (ACDM) cosmologies (e.g. Dave et al. (2010) ). However, because its density is relatively 
low, the detection remains a challenge (Dave et al. 2001; Dietrich et al. 2012; Werner et al. 2008). 

2. Giant radio galaxies (GRGs) 

Giant radio galaxies (GRGs) (linear size > 0.7 Mpc, h = 0.71) are a quite unique category, 
as they are the biggest objects in the sky with large radio structures. It is generally believed that 
their radio lobes are powered by beams of relativistic particles and magnetic fields from an active 
nucleus. Radio structures of many known GRGs are asymmetric over their mega-parsec extents 
(Fara et al. 2001; Machalski, Jamrozy & Zola 2001; Saripalli et al. 2005; Schoenmakers et al. 
2000; Strom et al. 2013) (see Figure 1). As the radio lobes expand, they interact with the ambient 
medium and can act as tracers of this barely detectable material. The large extent of GRGs provides 
a good opportunity to investigate the ambient gas. 

3. Radio galaxies and cosmology 

By 1950, several compact radio sources had been associated with distant galaxies (Bolton, 
Stanley & Slee 1949; Ryle, Smith & Elsmore 1950), and with the identification of the very strong 
radio source Cygnus A (Baade & Minkowski 1954; Smith 1951) with a z = 0.056 galaxy, the po¬ 
tential of radio sources for probing the Universe became apparent. This expectation was confirmed 
when 3C 295 was shown to be associated with a z = 0.461 galaxy (Minkowski 1960). Simulta¬ 
neously, radio source structure was found to be double lobed (Jennison & Das Gupta 1953), with 
sizes ranging from parsecs to hundreds of kilo-parsecs. 

Much effort went into identifying (and measuring the redshifts of) radio sources from reliable 
catalogues, especially 3C. High resolution radio maps (made by the One Mile Telescope, West- 
erbork and the VFA) tied down the radio source structures. The angular size measured could be 
translated into linear dimensions. Inverse Compton scattering of the radio emitting high-energy 
electrons becomes increasingly important at greater redshift. This effect could "snuff out" radio 
sources at high z (Miley 1968) before they could attain the dimensions found among the largest 
radio galaxies in the local Universe. 
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4C 73.08 
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Figure 1: WSRT map at 21 cm wavelength of the radio galaxy 4C 73.08. The brighter western lobe is more 
compact than its eastern counterpart. A similar effect is seen in NGC 6251 (Figure.2). 


NGC 6251 



Figure 2: WSRT map made at 92 cm wavelength of the radio galaxy NGC 6251. Like all of the GRGs, it 
is srongly polarized at long wavelengths. Note how unequal the distances between the ends of the two lobes 
and the host are. 


3 





GRGs as probes of ambient WHIM 


Bo Peng 


DA 240 



Figure 3: WSRT map made at 49 cm wavelength of the radio galaxy DA 240. Galaxies which are members 
of the DA 240 group have their positions indicated by symbols. 

4. Tracing the WHIM with GRGs 

As was clear from early WSRT observations of DA 240, 3C 236 and NGC 6251 (Figure 2), 
giant radio galaxies are strongly polarized at 608 MHz (Willis, Strom & Wilson 1974). The implied 
low ambient densities were at least qualitatively consistent with the large linear sizes of the radio 
components. A definitive study of five GRGs carried out by Mack et al. (1997) mapped their linear 
polarization at wavelengths from 92 to 2.8 cm. All of the sources were found to be polarized in 
many places at the longest wavelengths, indicating an internal rotation measure of ~ 1 rad nr~ 2 . 
The implied densities within the lobes depend upon the magnetic field strength. For equipartition 
field strengths of typically 1 fi G, the particle densities are of order 10 5 cm~ 3 . In a more general 
consideration of densities associated with GRGs, Mack et al. (1998) conclude that the IGM with 
which the radio components interact has a density of 1 - 4x 10 3 cm~ 3 . 

Several recent studies have explored the use of GRGs for probing the WHIM. In both the cases 
of Subrahmanyan et al. (2008), and Safouris et al. (2009), sensitive radio continuum images of two 
GRGs were used to relate their morphology to the ambient medium on the Mpc scale. They made 
spectroscopic observations to investigate galaxy distribution within a 2° field surrounding the radio 
galaxy hosts. Both studies showed clear correspondence between the giant synchrotron lobes and 
the surrounding galaxy distributions. Such a correspondence was also found in a series of papers 
(Chen et al. 2011b, 2012a,b). They found the existence of a galaxy group around each giant radio 
galaxy (NGC 6251, NGC 315 and 4C 73.08) through optical spectroscopic observations. However, 
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in the case of DA 240 (Chen et al. 2011a; Peng et al. 2004), the relatively symmetrical lobes 
coincide with a galaxy group. The group members evenly reside along the major axis of the lobes 
while the isolated host galaxy lies in the middle (Figure 3). The expected X-ray luminosities are 
much higher than the observed values, which suggests that the matter density within the lobes is 
quite low. A similar conclusion was obtained by Malarechi et al. (2013). In the investigation of a 
sample of 12 GRGs, they derived lobe energy densities from the radio observations via equipartition 
arguments. The inferred pressures in the lobes of the giant radio sources, which range from l.lx 
10~ 15 to 2.0 x 10 14 Pa, are lower than inferred from X-ray observations of dense filaments. This 
demonstrates the potential of using GRGs to study the WHIM gas within which their lobes evolve 
and interact. 

5. Prospects with SKA1 

A significant fraction of known giant radio galaxies were identified from the surveys of the 
WSRT and VLA: WENSS, NVSS and FIRST. The SKA1_SUR will be equipped with 96 dishes in 
Western Australia, including 36 dishes of ASKAP (latitude is about 27°S). At L band, it will have 
a high survey speed figure of merit (FoM) of about 156 times greater than the VLA, and 1.5 times 
better sensitivity as well as a 1.6 to 49 times higher angular resolution (Table 1). If it carries out a 
survey similar to the NVSS, with integration times of 10 minutes and a bandwidth of 500 MHz at 
L band, the noise level should achieve 9.1 pJy/beam. Since the GRGs usually have intermediate 
power between 10 25 and 10 26 W/Hz at 1.4 GHz, if the SNR is presumed (based on experience) 
to be 150 (Machalski, Jamrozy & Zola 2001; Schoenmakers et al. 2000), the luminosity distance 
will be 7825 Mpc (z ~ 1.2). Since the minimum elevation angle is about 15 degree (similar to 
the ASKAP), SKA1_SUR can observe the sky up to +48° (more than 80 percent of the whole 
sky). If the Galactic latitude is set to be \b\ > 12.5° to avoid the confusion of the Milky Way, 
the sky coverage will then be reduced to about 27,000 deg 2 . Adopting a number density of 10 7 
Mpc~ 3 (Saripalli et al. 2005), the SKA1_SUR is expected to detect about 12,000 GRGs in a total 
observing time of about 250 hours (formula listed in the Appendix). Even at an early phase (50% 
capability), SKA1_SUR will still be able to detect about 6,000 GRGs, which can be identified with 
counterparts obtained from optical surveys, such as SDSS, 2dF, and the planned LSST. The GRG 
sample will then increase by more than an order of magnitude. 

At the same time, the SKA1_MID will be equipped with 254 dishes in South Africa, including 
64 dishes from MeerKAT. Compared with the VLA, this will provide 6.2 times better sensitivity 
and 74 times higher survey speed. Moreover, the angular resolution will be 6.4 to 200 times better. 
After GRG candidates have been obtained from the SKA1_SUR, SKA1_MID can be used for 
a deeper investigation, especially for low surface brightness lobes, which is important for their 
use as cosmological probes. If the SKA2_MID could be equipped with PAFs, the survey speed 
with SKA2_MID would be greatly increased, and this will make the SKA2_MID quite a powerful 
instrument for discovering more of these faint, extended radio sources at higher redshift. 

Furthermore, since the frequency v of radio synchrotron emission is proportional to the square 
of the electrons’ energy E and the ambient magnetic field B (v ~ E 2 B), which is expected to be low 
and weak in the diffuse radio lobes of GRGs, detection of new radio structures in GRGs is likely 
from observations using the SKA1 low frequency array. There may even be discovery of new ones 
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Table 1: Summary of system sensitivity, survey speed and angular resolution of SKA1, VLA, WSRT at L 
band, and LOFAR at low frequency (extracted from Dewdney et al. (2013)). 


JVLA WSRT SKA1_MID SKA1_SUR LOFAR SKAl_LOW 


Asys/T S y S 

m 2 /K 

265 

124 

1630 

391 

61 

1000 

Survey Speed FoM 

deg 2 m 4 K 2 

1.76E4 

3.84E3 

1.3E6 

2.75E6 

5.21E4 

2.7E7 

Resolution 

arcsec 

1.4-44 

16 

0.22 

0.9 

5 

11 

FoV 

deg 2 

0.25 

0.25 

0.49 

18 

14 

27 

Bandwidth 

MHz 

1000 

160 

250 

500 

4 

250 


(like the discovery of GRG UGC 09555 with LOFAR). The low frequency array of SKA1 will be 
equipped with 250,000 antennas in Western Australia, with 16 times higher sensitivity and 520 
times quicker survey speed compared with LOFAR. 

Recently, Luo & Sadler (2010) modelled the evolution of FRI and II radio sources, finding that the 
latter evolve with decreasing monochromatic radio power, while the power of the former increases with age. 
The final stage for both sees rapidly declining emission as the radio source fades away, the lobes receiving 
little or no energy from the host nucleus or jets. As the lobes decline, their radio spectra will steepen, so a 
dying radio galaxy should have lobes with a steep radio spectrum and no jets or nuclear component. Finding 
candidates for this GRG phase will be well-suited to the SKA1 capability at low frequencies. 

With the highly improved system sensitivity, survey speed and angular resolution of SKA1, there should 
be more than one order of magnitude increase in the number of GRGs from about 140 at present (a small 
fraction of which have redshifts greater than 0.5 (Komverg & Pashchenko 2009)). With their ability to trace 
the ambient IGM, more GRGs will help illuminate the distribution and physical state of the WHIM and 
identify whether it can explain the missing baryon problem. And with more details on high-z GRGs, being 
able to compare them with low redshift sources will help to clarify how giant radio galaxies evolve, and 
especially how their environments affect their radio structures, and whether there are significant differences 
in group size, composition, morphology, etc., with cosmic epoch. 

At low redshifts, SKA1 with its excellent sensitivity and angular resolution will be able to investigate 
GRG groups in unprecedented detail. In several groups (like DA 240, Chen et al. (2011a); Peng et al. (2004); 
and 4C 73.08, Chen et al. (2012b); Strom et al. (2013)) we find that a large fraction of the companion 
galaxies are weak radio sources. The SKA1 sensitivity will enable us to explore the weak radio emissions 
of companion galaxies further. In addition, a number of the group galaxies appear to be spiral galaxies. 21 
cm HI observations of them are likely to turn up a number of detections. One can use the HI data to estimate 
the amount of mass present, and study the group kinematics. 


6. Conclusions 

Our goal with SKA1 is to carry out a large-scale survey to search for giant radio galaxies. The survey 
should be unbiased and complete down to a well-defined brightness or flux density, over a specified sky 
area. The GRGs found will be used to trace out the WHIM, while we also hope to detect candidate dying 
RGs, objects of low surface brightness with little evidence for compact features (jets, nuclear components), 
but with a steep radio spectrum. The high-z GRGs will be used to investigate how giant radio sources 
evolve. Some of the relatively nearby groups (like that surrounding DA 240) appear to have several spiral 
companions. We plan to observe them in the 21 cm HI line to study the motion and masses of GRG group 
members. 

Among the questions we hope to answer, we note the following: Can one use the GRGs to detect 
the WHIM? And if so, can we pinpoint the missing baryons? Can one use GRGs and their surrounding 
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companions to trace the distribution and physical state of the WHIM? Finally, are there enough high--; giant 
radio galaxies to investigate how radio galaxies evolve? 
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Appendix 

A. Formula used to estimate the 
number of GRGs for SKA_SUR 

Number of giant radio galaxies: 

N = p x Fco 

Number density: p = 10 -7 (Mpc -3 ) 

Comoving volume: (Mpc 3 ) 

Vco — vY^y/fcO 
Sky coverage: (deg 2 ) 


(A. 1) 


(A.2) 


Sobs = SNR x a (A.6) 


Sigma to noise: SNR >150 

Noise level with two polarizations: (Jy/beam) 


a = 


s/2Br 

System equivalent flux density: 

2^b Tsys 


Ssys — ' 


-Aeff 


(A.7) 


(A.8) 


Boltzman constant: k\, = 1.38 x 10 23 J/K 

For SKA1_SUR, sensitivity specification: 4^ = 

J sys 

391 m 2 K -1 ; bandwidth: B = 500 MHz 

Integrated time: 1 = 10 minutes with NVSS observ¬ 
ing model 


B. Estimate of the total observing time 


Asky = A RA x Asin(Dec) —Alx Asin(7>) (A.3) 

Comoving distance: (Mpc) 

Dl 

D co = -r-f- (A.4) 

1 +z 

Redshift z: (www.icosmos.co.uk) 


T A 


t = 


sky 


FoV 


(B.l) 


For SKA1_SUR, field of view: FoV = 18 (deg 2 ); 
integrated time: z = 10 minutes with NVSS observing 
model; sky coverage: A^y (equation A.3) 








